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We present a method for numerically obtaining the posi- 
tions, widths and wavefunctions of resonance states in a two 
dimensional billiard connected to a waveguide. For a rect- 
angular billiard, we study the dynamics of three resonance 
poles lying separated from the other ones. As a function of 
increasing coupling strength between the waveguide and the 
billiard two of the states become trapped while the width of 
the third one continues to increase for all coupling strengths. 
This behavior of the resonance poles is reflected in the time 
delay function which can be studied experimentally. 

In present-day high-resolution experimental studies, 
the properties of individual resonance states can be in- 
vestigated even when the level density is high. As an 
example, nuclear states have recently been identified and 
studied experimentally at very high excitation energy [Q . 
Collisional damping at such an energy is the same as at 
low excitation energy (ground-state domain). This ex- 
perimental result, being in contradiction to the standard 
statistical theory of nuclear reactions, can be justified 
by taking into account the interferences between reso- 
nance states arising from their interaction via the contin- 
uum 1^. In atoms, the coherent coupling of autoionizing 
states is studied experimentally e.g. in ||^. The authors 
point to the numerous possible applications of such in- 
vestigations. Thus, the interferences between individual 
resonance states play an important role and need to be 
considered in detail. 

Theoretically, such interferences have been studied in 
different fields of physics, see e.g. |4|-j^]. One important 
result is the phenomenon of resonance trapping, which 
arises from the interference of resonance states coupled 
to the same decay channel. As a function of a parameter 
controlling the coupling to the continuum, the widths 
of all states increase as long as the states are isolated. 



At a certain critical value of the parameter where the 
resonance states start to overlap each other, the widths 
bifurcate: the width of one of the states increases further 
while the width of the other ones decrease, see e.g. [|]. 
In other words: one of the states aligns with the channel 
and becomes short lived while the other ones decouple 
from the channel and become long-lived in spite of the 
strong coupling to the continuum. 

The phenomenon of resonance trapping is theoreti- 
cally well established but not proven directly up to now 
in experimental studies. In this letter we investigate 
the behavior of three neighboring resonances in a two- 
dimensional billiard connected to a single waveguide. As 
a function of the coupling between the resonator and the 
waveguide, we calculate both the position of the corre- 
sponding resonance poles and the Wigner-Smith time de- 
lay function. The time delay function, containing unique 
information on the interference between the resonances 
states, can be studied experimentally. 

The model used is as follows. We consider a two di- 
mensional billiard coupled to a waveguide (for instance 
a flat electromagnetic resonator). We have to solve the 
equation 



(-V^ + Ay)* = 



(1) 



where V is a potential barrier between the billiard and 
the attached lead. We use the Dirichlet boundary con- 
dition, ^E* = 0, on the border of the billiard and of the 
waveguide. The waveguide has a width W and the wave- 
function inside it has the asymptotic form 



^ = (e*'==" - R{E)e-''''') u{y) . 



(2) 



Here u{y) is the transversal mode in the waveguide, k is 
the wave number and R{E) is the reflection coefficient. 
It holds E — + (n/W)'^. We choose the energy of the 
incoming wave so that only the first transversal mode in 
the lead is open, i.e. 



u{y) = sin{—y) 



(3) 



Since \R{E)\ = 1, we write R{E) = exp{ie{E)) with 
Q{E) real. A measurable quantity derived from R is the 
Wigner-Smith time delay function. 
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Tm is the time the wave spends inside the billiard 

The energies and widths of the resonance states are 
given by the poles of the function R{E) analytically con- 
tinued into the lower complex plane. To find the poles we 
use the exterior complex scaling method Q . The general 
idea is to study the system after a scaling transformation 
applied to the x coordinate, see [0,||: x ^x = g{x). The 
function g is chosen as 
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X X > Xq 

9f{x) X < xq 



(5) 



with f{x) such that g{x) is three times smoothly dif- 
ferentiable and the inverse transformation g^^{x) exists. 
The attached waveguide extends from — oo parallel to the 
X-axis and we choose xq to be localized inside it. The re- 
lated transformation of the wavefunction reads 



^{x,y) 



1 



Using it the equation (Q) becomes 



(6) 
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^{i,y)^E^{i,y) (7) 



For a real parameter 6, this equation is fully equivalent 
to (|l|) since the transformation ^ is unitary. Moreover, 
the two equations are fully identical for x > xq. Since 
Xq lies inside waveguide the shape of the resonator is not 
changed by the transformation (||) which only rescales a 
part of the x-axis related to the waveguide. Moreover, 
since the waveguide is oriented parallel to the x-axis, the 
transformation does not change the boundary of the sys- 
tem. For 9 complex, (p) ceases to be unitary and the 
spectral properties of (l|) and are different. The 

continuous spectrum of (^) extends over ({tt/W)^ ,oo^, 

whereas the continuous spectrum of (Q) is rotated into 
the complex plane and is equal to 



U„=i:oo [{n^/Wf + e-^ (0,oo)} 



(8) 



This is a union of half-lines representing the continuous 
spectrum starting out from the real axis at every thresh- 
old energy (riTr/VF)^ with an angle —2 arg 9. The rotated 
continuous spectrum uncovers additive complex eigenval- 
ues of (0), the positions of which are independent of 9. 
These eigenvalues coincide with the resonance poles . 

In the following we study the time delay Tw and the 
resonance poles of a rectangular billiard of size Ax x 
Ay = 2 X 3.14 connected to a single waveguide with width 
W = 0.6. We choose y as a rectangular potential barrier 
with height 1 located at —0.3 < a; < 0. By changing 
the parameter A we can tune the coupling between the 
waveguide and the resonator. We calculate t^, by solving 
Eq. (|l|) with the Dirichlet boundary condition = 0, 
and the asymptotic boundary condition (||) imposed at 
X — —13. The resonance poles are found by the method 
of exterior complex scaling described above using xq = 
-2. 

In Fig. 1 we show the calculated time delay as a 
function of A and energy (panel a) as well as the depen- 
dence of the resonance poles on A (panel b). At large A 



(weak coupling to the waveguide) we see three isolated 
resonance states. As A decreases (the coupling to the 
waveguide increases) the lifetimes of all three states de- 
crease. The states attract each other in energy. As the 
resonances start to overlap, two of them become trapped 
while the third one gets short lived. At further decreased 
A the lifetimes of the two trapped resonance states in- 
crease. The motion of the poles is reflected in the time 
delay function. The lifetime of the short lived state is, 
at small A, so short that it practically disappears when 
plotting the time delay. The numerical errors in the dis- 
tances among the resonances are so small that the details 
of energy attraction and resonance trapping are stable. 

Using the method of complex scaling we can also 
study the wavefunctions of the resonance states (Gamow 
states). The interference between the resonances states 
leads to a mixing of their wave functions with respect to 
the eigenfunctions of the resonator (A oo). We illus- 
trate this phenomenon in Fig. 2, where the wavefunctions 
of the three resonance states are shown for A = 44, 23.5 
and 0, i.e. under the condition when (a) they are iso- 
lated, (b) very near to one another and (c) two of them 
are trapped while one is short lived. 

For A ^ 0, the amplitude of the wavefunction related 
to the resonance state which finally evolves into a short 
lived state is very small inside the resonator. This corre- 
sponds to a very small time delay, i.e. to a small probabil- 
ity of staying inside the billiard. The trapping of the two 
long lived resonance states occurs in two different ways. 
Both wavefunctions correspond to (almost) pure bound 
states of the billiard as long as A is large and become 
mixed when their distance in the complex energy plane 
to the other two states becomes small. At still smaller A 
one of the resonance states demixes and approaches again 
its pure shape. In accordance to that, its energy for small 
A approaches that for large A. The wavefunction of the 
other trapped resonance state remains however mixed. 
Accordingly, its position at small A differs from that for 
large A and its wavefunction is mixed with that of the 
short lived resonance state. 

The analysis presented in this letter shows that open 
billiards provide an excellent possibility to study the dy- 
namics of resonance poles in detail. A realization can 
be achieved by means of flat microwave resonators con- 
nected to a waveguide where the coupling strength to 
the channel can be varied by hand. Such systems allow 
therefore to investigate directly the formation of differ- 
ent time-scales (resonance trapping) by tracing the cor- 
responding time delay function. In particular, it is pos- 
sible to show the contraintuitive result that the lifetimes 
of certain resonance states increase with increasing cou- 
pling to the continuum. The results of such investigations 
will help in analysing high-resolution experimental data 
in various fields of physics. 
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Figure 1 

Contour and surface plot of log{Tw) (1-a) for a rectan- 
gular billiard (for parameters see text). The darker the 
plot is, the longer is the time delay. The motion of the 
corresponding resonance poles with A (l.b). The posi- 
tions of the resonance poles for A = 44 are denoted by 
squares and for A = 23.5 and A = with large dots. 

Figure 2 

The wavefunctions of three resonance states for A = 44 
(case a)), A = 23.5 (case b)) and A = (case c)). The 
energies of the states at A = 44 are 38.6 (state N- 1), 
38.8 (N^ 2) and 39.8 (N^ 3). The location of the poten- 
tial barrier is marked by a rectangular in the waveguide. 
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